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Abstract 

Multiple production of Higgs particles is essential to study Higgs self-couplings 
at future high-energy colliders. In this paper we calculated the resonance con- 
tributions to the production of three lightest neutral supersymmetric Higgs 
bosons in gluon fusion at LHC. The cross sections due to trilinear Higgs cou- 
plings is sizeable but the measurement of the quartic coupling XhhhHih) seems 
to be impossible. 



1 Introduction 



One of the basic open questions of particle physics is the nature of electroweak symmetry 
breaking. Beyond the discovery of the Higgs boson(s), however, reconstructing the Higgs 
potential will be necessary, and that requires the experimental study of the self-couplings 
of the Higgs bosons. 

Pair production of neutral Higgs particles in gluon fusion sensitive to the trilinear 
couplings of Higgs bosons was studied in the SM [1] and the Minimal Supersymmetric 
Standard Model (MSSM) [2-5]. QCD corrections were included in the limiting case of a 
heavy top mass |J. 

In the SM the trilinear and the quartic couplings of the physical Higgs particle are 
fixed by the Higgs mass. There is no resonance in the pair production of SM Higgses, 
and the cross section is only ~20-50 fb in the intermediate mass range [3]. No resonance 
effect is present in the production of three Standard Model Higgs bosons via one Higgs 
in gluon fusion and it is estimated to be under the discovery limit at LHC. 

In the MSSM there are five Higgs bosons (h, H, A, H^) and many trilinear and quartic 
couplings among them. Two parameters describe the Higgs sector of the MSSM: Ma, the 
mass of the CP-odd Higgs boson A, and tan/?, the ratio of the two vacuum expectation 
values. For a wide range of Ma, in pp collision only the cross section of the hh production 
is large [2,3]. Possible other processes, such as the gauge boson fusion and Higgs-strahlung 
off W bosons or heavy quarks (top) provide less event number. 

In order to provide further tests of trilinear and quartic Higgs self-couplings, in this 
paper we calculate the resonance enhanced production of three lightest supersymmetric 
Higgs bosons (h) in gluon fusion 

PP~> 99 -> hhh. (1) 

Generally four different graphs contribute to the production of three Higgses. Gluons 
couple to triangle, box or pentagon quark loops emitting 1,2 or 3 Higgs bosons, as is seen 
in Fig.l. Squark loops are neglected. We have omitted the pentagon graph producing a 
flat continuum background for the resonance production and we will give the complete 
cross section valid in a wider range in a subsequent publication |7j]. We learn that the 
quartic coupling is not accessible by LHC experiment. The resonance contribution to (1), 
however, sizeable at LHC For instance at tan (3 = 3 it yields a cross section of about 300 
fb. This is about 1/10 times the corresponding hh-ca.se [3]. 

The masses, widths and the couplings were calculated using the complete one-loop 
and the leading two-loop radiative corrections from || . The relevant Higgs self-couplings 
are the following 



3e cos 3 a 

Ahhh = 3 cos(2a) sm(p + a) + — j , 

M z sin p 

36 sin cy. cos^ ex 

^Hhh = 2 sin(2a) sin(/5 + a) — cos(2a) cos(/5 + a) + 

IVl 7 



sin (3 



2 



. . . N . . . . . . . 3e sin 2 a cos a 

Affif/i = —2 sin(2a) cos(p + a) — cos(2a:) sm(p + a) + —f " ~ 



M| sin /3 



Anfehh = 3 sin(2a) cos(2a) + 



e cos a sin 2 a 



8M| sin 2 /3 



sin 4 a 



W =3cos( 2a ) + i ^-^ (2) 

The trilinear and quartic couplings are normalized to A3 = [v / 2G , i?] 1//2 M| and A4 = 
\[2GpM\, respectively. The couplings depend on (3 and the mixing angle a of the CP- 
even Higgs sector 

tan 2a = — = -j, f- — tan 2(3. (3) 
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The common squark mass is fixed to Mg= 1 TeV. The quartic couplings depend on M4 
and are shown for two values of tan/3 = 3 and 30 in Fig. 2. 

The paper is organized as follows. In the next section we give the cross section of the 
resonance production of three h's in gluon fusion and in section 3 we show the numerical 
results. 



2 The cross section 

Three different channels contribute to the resonance production of three lightest supersym- 
metric Higgs particles (h) in gluon fusion, (a) One virtual Higgs boson (h,H) is produced 
by the heavy quark triangle and it decays into 3 h's via the quartic coupling \Hhhh or 
\hhhh, Fig- la. (b) One virtual Higgs boson decays into 3 h's in two steps testing trilinear 
couplings, gg — > H, h — > (H, h)h — > hhh, Fig. 2b. (c) Heavy quark box diagram couples 
to two Higgses, one of them decays subsequently into two h's, gg — > (if, h)h — > hhh. Fig. 
2c. We have omitted the pentagon graphs contributing only to the continuum production. 
We believe it increases the production at large tan/3 due to the large hbb-coupling, but 
does not change our conclusion about the measureability of the quartic Higgs-couplings. 

There are two different helicity amplitudes contributing to the total cross section: the 
total spin of the two gluons along the collision axis can be S z — (contribution F) or 
S z = 2 (contribution G). The triangle graphs give only contributions F 3 (S z = 0) as they 
involve a single spin-0 Higgs intermediate state. Box graphs give both contributions F4 
(S z = 0) and G4 (S z = 2). F 4 and G 4 are Lorentz and gauge invariant decompositions of 
the box amplitude. The tensor basis is the following 
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where pi,P2 are the momenta of the incoming gluons and p^ is one of the momenta of the 
outgoing Higgs bosons {ps, P4, P5) ■ Here p\ = 2 ^ Pl ^^ P3 ^ —p\ is the transverse momentum 
of the third particle. Tensors A^ v and are orthogonal and normalized to 2 [2]. 
The M-matrix of the process is 
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The spin and color averaged parton level cross section of the process gg — > hhh is 
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Here we used the Chew-Low parametrization of the three particle phase space ||10|| , a 
trivial angle integration was carried out, §ik = (pi + Pk) 2 , and the usual lambda function 
is X(x,y,z) = (x — (y/y + i/i) 2 )(a; — (^/y — V^) 2 )- dVt^ 1 is the differential solid angle 
in the center of mass system of particles 4 and 5. 

F receives contributions from the triangle and box graphs, while only box diagrams 
give contributions S z — 2 to G 
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Here we separated the functions F3, F4, G4 responsible for the heavy quark triangle and 
boxes, C 3 , C 4 contain the coupling constants and the propagators. The amplitude has 
to be summed over all quark flavours. However, the four light quarks can be neglected 
having very small Higgs couplings. Supersymmetric particles are assumed to be too 
heavy to participate in the triangle or box loop. The second sum in the box contribution 
corresponds to the outgoing Higgs particle that couples directly to the quark loop. 

We calculated F 3 , F4, G4 and found them in agreement with the results of 0, For 
instance, 
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F4, G4 have long analytic expressions which can be found in ref 0. The box functions 
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In case of a large quark mass {rrvq ^> s ~ Mf lH ) we get as in [2] 

F 3 = 2 - + 0(s/ml), if = -| 0(s/ml), Gf = O (i/m 2 Q ) (11) 
In the limit of light quark masses (mg <C s ~ ^hu) a ^ the form factors vanish to 

o( m ys). 

The generalized triangle couplings are the following 

C 3 = C% + + Cf + C% H + Cf h + Cf H 

M 2 

C *' = XhhhHa s-Ml a AM Ha T Ha 9 Q (12) 

t> - *H a H b h- T72 1 - n/r ^ *H b hh- 772 1 - n/r 9q , tta,b ~ fl tt. 

s ~ M L + lM nSn a s - Mfj b + iM Hb T Hb 
The box couplings are 

4 — t - / 4 "T ^4 

r (h,H){i) _ x „fc„ fc / fl 

SfeZ — ■ /M ? l /H + lM h/H*- h/H 

i, k, I are cyclic. gA denotes the Higgs-quark couplings normalized to the SM Yukawa 
coupling [V2G F ] 1/2 m Q , 

= cos a/ sin/?, g% = — sin a/ cos/3, = sin a/ sin /5, = cos a/ cos/5. (14) 

3 Results 

Now, we calculate the production cross section of three lightest neutral supersymmetric 
Higgs bosons pp —>■ hhh + X via gluon fusion at the LHC energy of y/s = 14 TeV in the 
above picture by folding the parton level fusion cross section with the gluon luminosity. 



We used the GRV structure functions [11] for the gluon luminosity at Q = s. The 



numerical integration was made by the VEGAS package ||12|| . The total cross section vs. 
Ma is shown in Fig. 3 for two represantative values of tan/5 = 3 and tan/3 = 30. In the 
plotted range of Ma the heavy CP-even Higgs boson (H) is nearly degenerate in mass 
with A while the mass of the lightest neutral Higgs h is quickly approaches approximately 
104 GeV from below. 

For tan/5 = 3 in the range of Ma = 200-350 GeV the cross section is enhanced by 
resonance effect and large. The cross section goes up to nearly 300 fb and slowly decreases 
to about 150 fb giving a large number of events at the expected luminosity / Ldt ~ 
100 fb^ 1 /year. The main contribution comes from the box diagram where the heavy 
CP-even Higgs boson H couples to the quark loop and decays into two bosons h. Mh 
slowly increases around 100 GeV when Mh — Ma reaches 200 GeV the H propagator 
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enters the resonance region. At M H ~ 2m t there is a threshold effect due to the fall-off 
of the branching ratio BR(H — > hh) and partly because on-shell top quark pairs can be 
produced in the quark loop. For large values of Ma the cross section reaches a continuum 
value of a few fb. The dip at small Ma, similarly to the case of tan/5 = 30, is the result 
of the zero in the trilinear couplings X(H/h)hh- 

For tan/3 = 30 the main contribution comes from the box diagram where an h couples 
to the quark loop, propagates and decays into two h's via Xhhh- The Hhh coupling 
is small compared to the case of tan/5 = 3 while the hhh coupling is larger giving a 
sizeable continuum. There is no observable resonance effect and the cross section is 
nearly constant, 6 fb, versus M A after the coupling constant changed sign. The K factor 
is expected to increase the cross section as in other cases [6,3]. 

The contributions of the diagrams containing a heavy quark triangle are also shown in 
Fig.3 at tan/3 = 3. The lowest curve corresponds to the diagram with quartic couplings, the 
one in the middle sums up all the contributions of the graphs involving a quark triangle. 
There is a clear resonance effect in both curves when M H reaches 3M h at ~310 GeV. Here 
the first propagator of gg — > H — > Hh — > hhh becomes resonant too. The first rise in the 
middle curve is the result of the resonance in the second propagator of Fig. 2b, similarly 
to the box. The triangle loops yield, however, only a small fraction of the cross section 
even at their peak at 300-350 GeV. 

The suppressed contribution of the quartic coupling implies that in pp — > hhh + X 
the measurement of the quartic coupling is not possible. Not only the box contribution 
is larger by two orders of magnitude but also the other triangle diagrams are 10 times 
larger. For tan/3 = 30 the case is even worse. 

In conclusion, in this paper we have calculated the resonance contribution to the cross 
section of 3h-production at LHC via gluon fusion. The main contribution comes from the 
trilinear Higgs coupling while the quartic ones turn out to be negligible. In the resonance 
region at small tan/3 the ratio of 2h and 3h production is about 10. 

This work was partially supported by Hungarian Science Foundation Grants under 
Contract Nos. OTKA T029803 and F022998. 
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FIGURE CAPTIONS 



Fig.l: Generic diagrams contributing to the production of three CP-even MSSM Higgs 
bosons in gluon-gluon collisions, gg — > hhh: a) triangle with quartic coupling, b)triangle 
with trilinear couplings c) box , d) pentagon contributions. 

Fig. 2: Quartic Higgs couplings in the MSSM as functions of the pseudoscalar Higgs mass 
Ma for two representative values of tan/5 = 3 and 30. 

Fig. 3: Total cross sections for production of three lightest CP-even MSSM Higgs bosons 
in gluon-gluon collisions at the LHC for tan/3 = 3 and 30 (upper two curves). The upper 
axis presents the scalar Higgs masses Mh for tan/5 = 3 corresponding to the pseudoscalar 
masses M4. The lower two dashed curves represent the resonance contributions of the 
triangle graphs of Fig. 2a and 2b for tan/5 = 3. 
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